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Abstract 

The high-pressure behaviour of a-PbF2 has been investigated by angular-dispersive 
synchrotron X-ray powder diffraction up to 6.55(8) GPa. The fit of a 3 rrf -order Birch- 
Murnaghan equation-of-state gave the volume at zero pressure Vo=194.14(4) A 3 
and the bulk modulus bo=47.0(6) GPa with the pressure derivative b'=7.9(4). The 
continuous-symmetry-measure approach has been used for the quantification of the 
distortion of the coordination polyhedron in a-PbF2 revealing an increasing distor- 
tion with pressure. 
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1 Introduction 



Lead fluoride can be found in two polymorphs at ambient conditions: /3-PbF 2 
in the fluorite type structure [1] and a-PbF 2 in the orthorhombic cotunnite 
type structure (fig. 1) with space group Pnma and a=6.440 A, 6=3.899 A, and 
c=7.651 A [2]. The structure is built up by corrugated layers of tricapped trig- 
onal PbFg prisms along a sharing an edge parallel to the b direction. All atoms 
occupy the Wyckoff position 4c (x, 1/4, z) with x=0.2527(13), z=0.1042(7) 
for Pb, x=0.8623(21), z= 0.0631(15) for F(l) and x=0.4662(20), z=0.8457(17) 
for F(2), respectively [2]. 




Fig. 1. Schematic drawing of the a-PbF2 crystal structure. The black spheres repre- 
sent Pb 2+ and the gray ones F~ . The first coordination shell around Pb 2+ is shown 
where at F(l), (5), (8), (9) and F(2), (3), (4), (6) and (7) occupy symmetrically 
equivalent positions. 

It is not clear which polymorph is thermodynamically stable at ambient con- 
ditions. In recent years both phases have intensively been studied due to their 
potential industrial applications [3,4,5,6,7,8,9,10]. /3-PbF 2 shows a high supe- 
rionic conductivity with a low superionic transition temperature suitable for 
solid state batteries and sensors [5]. The a-phase is a strong scintillator useful 
in detector technology [7] . The high-pressure behaviour of PbF 2 has been stud- 
ied recently by various experimental and theoretical techniques. Lorenzana et 
al [8] investigated the high-pressure properties of PbF 2 by Raman spectroscopy 
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in conjunction with a first principles, full-potential linear muffin-tin orbital to- 
tal energy method in the pressure range from ambient pressure up to 30 GPa. 
The authors detected a phase transition at approximately 0.5 GPa from the cu- 
bic /9-PbF 2 to the orthorhombic a-PbF 2 . At 14.7(5) GPa a second phase tran- 
sition to 7-PbF 2 with a yet unknown structure has been observed. Considering 
previously investigated high-pressure phases of isotype ionic AX 2 compounds 
(e.g. SnCl 2 , PbCl 2 , BaCl 2 and BaF 2 [11,12,13]), Lorenzana et al proposed a 
hexagonal Ni 2 In type structure and a monoclinic structure with a doubled a 
axis as possible structures for 7-PbF 2 . Their calculations revealed that above 
16.4 GPa, PbF 2 in the Ni 2 In type structure is more stable than a-PbF 2 . How- 
ever, 7-PbF 2 in a hexagonal symmetry could not explain the observed Raman 
spectra and consequently, the authors proposed the monoclinic structure for 
the high-pressure 7 polymorph. High-pressure and high-temperature neutron 
diffraction experiments have been performed from ambient conditions up to 
pressures of 1 GPa and temperatures of 966 K to investigate the phase diagram 
and to determine the nature of the disorder of the fluorine ions at high tem- 
perature and its pressure dependence within the a and (3 polymorphs of PbF 2 
[9]. The phase stability of PbF 2 under high-pressure has also been studied by 
ab initio calculations using the perturbed ion method [10]. The calculations 
have been performed for the a and j3 polymorphs as well as for the two pro- 
posed 7-PbF 2 structures. The results of the calculations support the existence 
of 7-PbF 2 in the hexagonal Ni 2 In type structure above 20.2 GPa. Experimen- 
tally determined high-pressure data are mainly for cubic /9-PbF 2 [8,9]. Hence, 
the goal of the work presented here was to measure the structural changes 
induced by high-pressures in order to elucidate the compression mechanism 
of a-PbF 2 and to extend the limited pressure range of the study of Hull and 
Keen [9] for a precise determination of the equation-of-state. Here we present 
the results of a high-pressure synchrotron powder diffraction study of a-PbF 2 
up to 6.55 GPa. 



2 Experimental 



The X-ray diffraction experiments were performed using a commercial (Merck, 
omnipur) polycrystalline sample of PbF 2 . By conventional X-ray powder diffrac- 
tion it was found that the sample is pure a-PbF 2 . High-pressure powder 
diffraction experiments up to 6.55(8) GPa were performed at the wiggler beam- 
line 5.2R at the ELETTRA synchrotron in Trieste, Italy. The beam-line optic 
consists of a Si(lll) double-crystal monochromator in non-dispersive config- 
uration followed by a three-segment platinum-coated toroidal focusing mirror 
with a horizontal acceptance of 2.8 mrad. The incident beam was collimated 
to 80 /im diameter. High-pressure powder diffraction patterns were collected 
at the wavelength A=0.6888 A using a marresearch (mar345) image plate de- 
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tector with a pixel-size of 100 xlOO /im. The exposure time per image was 
about 20 minutes. Pressure was applied using a Diacell DXR-6 diamond anvil 
cell mounted on a rrz-positioning-stage for the adjustment in the incident 
beam. The sample was placed in the hole (0 =150 /im) of an Inconel gas- 
ket preindented to 100 /im. In order to ensure hydrostatic conditions a 16:3:1 
methanol-ethanol- water mixture was used as pressure transmitting fluid [15]. 
We used the ruby fluorescence method for the pressure determination applying 
the Piermarini pressure scale [16]. The sample to detector distance was deter- 
mined by measuring the sample in the uncompressed cell. Geometry param- 
eters for the radial integration of the two-dimensional data were determined 
using FIT2D [17]. For the transformation into standard one-dimensional pow- 
der patterns the software TW020NE [18] was used, which allows for counting 
statistics of multiply measured data points in order to provide proper error 
estimates of the intensities and hence, a correct weighting scheme in sub- 
sequent least-squares refinements [19]. The one-dimensional powder patterns 
were corrected for attenuation by the pressure cell. Lattice parameters were 
obtained from profile matching and the structure was refined at five different 
pressures by the Rietveld method employing the program FULLPR0F [20] . The 
background was described by a linear interpolation between selected points 
and the peak profile was modelled using a pseudo-Voigt function [21]. The 
structural parameters from Boldrini and Loopstra [2] were taken as the ini- 
tial values for the refinements. The isotropic displacement parameters (Uj SO ) 
for lead and the two fluorine atoms were determined from the measurement 
at ambient pressure and assumed to be constant in the observed pressure 
range. The standard deviations of the refined parameters were scaled with the 
Berar-factor [22]. The pressure- volume data were fitted by a 3 rd -order Birch- 
Murnaghan equation-of-state [25]. Since pressure and volume are subject to 
experimental errors, which have to be weighted correctly in the fit procedure, 
the conventional least-squares fit is not appropriate. Consequently, the method 
of orthogonal distance regression (ODR) was used and the equation-of-state 
implemented in the program code 0DRPACK [23,24] to minimise 

Ew:-K)) 2 +K.^- y , c )) 2 . (!) 

i 

with weights w pi = C\loio and wyi = c 2 /o"yo. The parameters c 12 are ad- 
justable parameters ensuring that the weighted experimental errors w are in 
the same order of magnitude for both variables [24] . , The difference between 
the observed and the calculated values for pj and Vj are at least one order 
of magnitude smaller than the experimental error bars, when the experimen- 
tal error bars are used as weights in the fit-procedure. This indicates that 
the relative precision of the observed data points, which influences the fit, is 
higher than that represented by the error bars. Taking the differences Ap and 
AV between the observations and the equation-of-state as weighted errors in 
the subsequent fitting cycle did not affect the values of the fitted parameters 
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themselves but their error estimates. 



3 Results and Discussion 



As an example the observed and calculated diffraction patterns resulting from 
the Rietveld refinement of the 3.8(1) GPa data are presented in Figure 2. The 
pressure dependence of the normalised lattice parameters a/ao, b/bo, c/cq and 
the unit cell volume are shown in Figures 3 and 4, respectively. In Table 1 the 
lattice parameters and the residuals resulting from the Rietveld refinement 
at five selected pressures are given. The refined atomic positions and the re- 
sulting bond distances are given in Tables 2 and 3, respectively. The atomic 
displacement parameters (Uj SO ) were refined to 0.027(3) A 2 and 0.029(2) A 2 
for lead and the two fluorine ions, respectively 




Fig. 2. Observed (circles) and calculated (solid line) diffraction pattern for a-PbF2 at 
3.8(1) GPa. The difference and the tick marks for the calculated reflection positions 
are plotted at the bottom of the figure. 

The compression of a-PbF 2 is extremely anisotropic. The compression parallel 
to the a axis is 1.5 times larger than that parallel to the b axis and twice as large 
as along the c direction. The volume at zero pressure Vo=194.14(4) A 3 , the 
bulk modulus b =47.0(6) GPa and its pressure derivative b'=<9b/<9p=7.9(4) 
were determined from the equation-of-state fit. 
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Fig. 3. Evolution of the normalised lattice parameter with pressure. The compres- 
sion in the three lattice directions a (■), b (•) and c (A) is anisotropic. The ex- 
perimental error bars for the pressure and the lattice parameters correspond to 
the size of the symbols. The lines are curve fittings of 3 rd -order Birch-Murnaghan 
equation-of-states to the data. 

These results compare to Vo= 198.48 A 3 , b =40 GPa and b'=4.8 derived from 
quantum mechanical calculations at the density functional theory level [8] and 
with bo=57.87 GPa and b'=4.98 derived from Hartree-Fock calculations [10] 
within the accuracy usually obtained in such types of calculations [26,27]. The 
bulk modulus for a-PbF 2 of 117(4) GPa reported by Hull and Keen [9] differs 
significantly from the value derived in this work. However, this difference can 
be attributed to the limited pressure range and the linear fitting of the volume 
data for the determination of the bulk modulus in the work by Hull and Keen 
[9]. The values for the b and b' obtained in this work are in the same order 
of magnitude as found for other halides crystallising in the cotunnite type 
structure, e.g. PbCl 2 with b =34(l) GPa and b'=7.4(6) [13] and CaCl 2 with 
ab =51 GPa [14]. 

The changes of the atomic positional parameters in the investigated pressure 
range are small. The Pb-F distances in the tricapped trigonal prism range from 
2.38(2) to 3.12(3) A at ambient conditions (Tab. 3). The mean of the Pb-F 
distance in the coordination polyhedron (2.66(2) A) is close to the sum of the 
ionic radii [28] for Pb 2+ (1.35 A) in ninefold- and F - (1.33 A) in fourfold- 
coordination. Under a pressure of 6.55(8) GPa the average Pb-F distance is 
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Fig. 4. Pressure dependence of the unit cell volume of a-PbF2- The solid line results 
from the fit of a 3 r<i -order Birch-Murnaghan equation-of-state, dashed lines represent 
the 95 % confidence region of the fitted equation-of-state. The experimental error 
bars correspond to the size of the symbols. The insert shows the one-sigma error 
ellipsoid of bo and b' for the fit. 

reduced to 2.59(2) A. The individual Pb-F bonds show a different compression 
behaviour. The following notation of the fluorine ions is according to fig. 1. 
The bonds in the a6-plane decrease by 3.64 %, 7.34 % and 0.8 % for Pb-F(l), 
Pb-F(3) and Pb-F(8), respectively . The Pb-F(2) bond is increased by 4.62 %, 
Pb-F bonds oriented along the c direction are less compressible. The Pb-F(5) 
and the Pb-F(6) bonds are shortened by 1.8 % and 2.88 % at 6.55(8) GPa. 

The structural changes in a compound under high-pressure can be described by 
considering the distortion of the respective coordination polyhedron. Recently, 
two different approaches have been developed [29,30] to quantify the distortion 
of a polyhedron by comparison of an ideal to the real one. 

The capped square anti-prism and the tricapped trigonal prism (TTP) are 
ideal polyhedra for a nine-fold coordination [31]. The TTP is found in a- 
PbF 2 . The regular TTP, shown in fig. 5, is characterised as follows: all nine 
ligands are located on the surface of a sphere with radius r s , and the polyhe- 
dra edge lengths are AB=AC=2r s /\/3 and AD=2r s /3^5 [31]. The distortion 
measure of Mackovicky and Balic-Zunic [29] indicates the volume difference 
between real and ideal coordination polyhedra. Therefore the ligands of the 
real and ideal polyhedra have to be located on a common circumsphere with 
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Table 1 

Residuals, lattice and profile parameters from Rietveld refinement for selected pres- 
sures. 



p [GPa] 


0.00(4) 


1.62(5) 


3.80(9) 


5.44(9) 


6.55(8) 




11.1 


15.0 


14.1 


11.3 


13.6 


x 2 


1.93 


3.44 


3.87 


2.53 


2.16 


a [A] 


6.4567(9) 


6.372(2) 


6.278(1) 


6.213(1) 


6.165(1) 


b[A] 


3.9071(5) 


3.8680(9) 


3.8294(8) 


3.8088(8) 


3.7963(8) 


a [A] 


7.666(1) 


7.603(2) 


7.549(2) 


7.520(1) 


7.503(2) 


V 


0.13(4) 


0.07(3) 


0.25(4) 


0.14(3) 


0.09(4) 


u 


-0.19(9) 


-0.17(5) 


-0.85(2) 


-0.53(1) 


-0.47(3) 


V 


0.14(2) 


0.16(7) 


0.39(3) 


0.29(3) 


0.27(4) 


w 


0.0034(3) 


0.006(2) 


0.017(2) 


0.010(2) 


0.009(5) 



Table 2 

Atomic positions resulting from Rietveld refinements. 







Pb 


F(l) 


F(2) 




X 


z 


X 


z 


X 


z 


0.00(4) [GPa] 


0.254(3) 


0.1067(5) 


0.891(7) 


0.051(7) 


0.489(8) 


0.850(6) 


1.62(5) [GPa] 


0.253(4) 


0.1111(7) 


0.880(9) 


0.03(2) 


0.49(1) 


0.850(9) 


3.80(9) [GPa] 


0.251(5) 


0.1139(8) 


0.89(1) 


0.02(1) 


0.50(1) 


0.854(9) 


5.44(9) [GPa] 


0.253(4) 


0.1145(6) 


0.888(9) 


0.041(7) 


0.477(9) 


0.853(8) 


6.55(8) [GPa] 


0.254(4) 


0.1163(7) 


0.884(9) 


0.037(9) 


0.51(1) 


0.843(9) 



radius r s . However, this measure is not applicable to polyhedra with a large 
variation in the bond distances as found in a-PbF 2 . Hence, the distortion 
measure based on the continuous symmetry measure (CSM) [30] was used. In 
this method a regular tricapped trigonal prism is embedded in the observed 
one. As an additional free parameter the radius around the three cap atoms 
r e was used to consider non-spherical deformation. In order to find the best 
conformation, an evolution strategy [32] was used to minimise the euclidean 
norm (instead of the standard least-squares method employed by Pinsky and 
Avnir [30], which requires the computation of the two-dimensional derivative 
with respect to r s and r e ). For various pressures the results of the embed- 
ding are given in table 4. At ambient pressure the TTP is slightly distorted 
(CSM=1.378) compared to the ideal value (CSM=1). The nine fluorine ions 
are not arranged on the surface of a sphere (j q =l) but on an ellipsoid of rev- 
olution (r g =1.062). Remarkably, the lead ion is not located in the conjoined 
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Table 3 

Bond distances in the coordination shell around a lead atom (see fig. 1. The values 
are given in A. 



p [GPa] 


0.00(4) 


1.62(5) 


3.80(9) 


5.44(9) 


6.55(8) 


Pb-F(l,9) 


2.47(2) 


2.40(1) 


2.35(2) 


2.39(2) 


2.38(2) 


Pb-F(2) 


2.48(2) 


2.53(2) 


2.51(3) 


2.41(3) 


2.60(3) 


Pb-F(3,4) 


2.59(3) 


2.51(3) 


2.49(3) 


2.55(2) 


2.40(2) 


Pb-F(5) 


2.77(2) 


2.78(3) 


2.88(2) 


2.72(1) 


2.72(3) 


Pb-F(6,7) 


3.12(3) 


3.10(2) 


3.06(3) 


2.98(2) 


3.03(3) 


Pb-F(8) 


2.38(2) 


2.44(2) 


2.36(2) 


2.36(1) 


2.36(2) 



Table 4 



Evolution of the tricapped trigonal prism coordination polyhedra in a-PbF 2 with 
pressure. Tabulated are: The mean fluorine distance d m F from the centroid, the 
scattering of the fluorine distances around the mean a, the distance Pb exc between 
the lead ion and the centroid of the nine fluorine ions, the minimal distance d and 
the mean distance per coordinate d/coord between observed and embedded ideal 
polyhedra. The distortion is given in the continuous symmetry measure CSM [30] 
and the quotient r g =r s /r e of the radii around the circumscribed sphere (r s ) and 



caps (r e ). 


p [GPa] 


dmF [A] 


a [A] 


Pbexc [A] 


d[A] 


d/coord [A] 


CSM 




0.00(4) 


2.6574 


0.557 


0.3049 


0.8803 


0.033 


1.378 


1.062 


1.62(5) 


2.6324 


0.564 


0.3002 


0.9350 


0.035 


1.492 


1.024 


3.80(9) 


2.6098 


0.598 


0.3023 


1.1464 


0.042 


1.859 


1.003 


5.44(9) 


2.5857 


0.477 


0.2458 


0.8787 


0.033 


1.455 


1.052 


6.55(8) 


2.5767 


0.627 


0.2849 


0.9667 


0.036 


1.607 


1.006 



centroid of the fluorine ligands. The variance of the fluorine bond distances 
increase with pressure, reflecting the difference in the compression of the in- 
dividual Pb-F bonds. The lead ion shifts towards the TTP-centroid, and the 
arrangement of the fluorine ligands becomes more spheric with pressure. How- 
ever, the minimal distance between the observed and the embedded fluorine 
ions increases. Compared to TTP the distortion is larger and the spheroidal 
arrangement of the fluorine ligands extends, both effects can be attributed to 
increasing angular distortions. 

Raman spectroscopy measurement above 14.7 GPa in PbF 2 indicate a struc- 
tural phase transition to a post-cotunnite structure with monoclinic symmetry 
[8] as observed in the isotype SnCl 2 and PbCl 2 compounds [13]. The extrap- 
olation of the CSM to higher pressures indicate a increasing distortion of the 
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Fig. 5. The regular tricapped trigonal prism with all ligands located on a sphere with 
radius r s . The polyhedral edge lengths are AB=AC=2r s /\/3 and AD=2r s /3v / 5- 

TTP polyhedron. However, from the present CSM data a saturation of the 
CSM may not be excluded. Such a behaviour would suggest a ridged-unit like 
behaviour of the TTP as it can be found e.g. for the tetrahedral building 
units in silicates. Since the structure of the post-cotunnite phase is not known 
for PbF 2 , a further examination of the high-pressure behaviour seems to be 
justified. 
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